Impact of inertia on biased Brownian transport in confined geometries 
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We consider the impact of inertia on biased Brownian motion of point particles in a two-dimensional channel 
with sinusoidally varying width. If the time scales of the problem separate, the adiabatic elimination of the 
transverse degrees of freedom leads to an effective description for the motion along the channel given by 
the potential of mean force. The possibility of such description is intimately connected with equipartition. 
Numerical simulations show that in the presence of external bias the equipartition may break down leading 
to non-monotonic dependence of mobility on external force and several other interesting effects. 



Particle transport in channels attracted recently much 
attention due to its importance in zeolites^, biological^ 
and designed nanopores^, and other situations. This 
activity was stimulated by the interest to shape and 
size selective catalysis^, particles' separation^ and to 
polymer translocation^. The progress in experiments 
boosted theoretical activities^. Works by Jacobs and 
Zwanzig& lead to the so-called Fick- Jacobs (FJ) ap- 
proach, in which the elimination of transversal degrees of 
freedom (assuming fast equilibration in transversal direc- 
tions) results in an effective description for the longitudi- 
nal coordinate evolving in the potential of mean force. 
The approach found its application for biased trans- 
port through periodic planar—^ and three-dimensional^ 
channels with smoothly varying cross-section. Other 
approaches were applied to transport through cylin- 
drical septate channels^, tubes formed by spherical 
compartments^ and channels with abrupt changes in 
cross-section^. 

The applicability of the FJ approach depends on the 
existence of a hierarchy of relaxation times governed by 
the geometry of the channel and by friction. This hierar- 
chy guarantees the time scale separation and the equipar- 
tition of energy both being necessary conditions for ap- 
plying the method. The core physical assumption behind 
the FJ approach is that the the dynamics of particles in 
a fluid inside the channel is the overdamped Langevin 




Figure 1. (Color online) Sketch of a segment of a sinusoidally- 
shaped 2D channel and an example for a particle trajectory 
for 7 = 1 and / = 1. The constant force / is pointing in longi- 
tudinal (x-) direction. Shown are also the potentials U±(n±) 
of the particle-wall interaction. 



the Langevin equation 



-7v-V q f/(x,y) + V27^(*) 



(1) 



oneii. In this work we study the impact of the friction 
coefficient 7 and of external bias / on the transport in the 
channel and concentrate on its influence on equipartition. 

We consider Brownian particles with mass to in a 2D 
channel with top and bottom boundaries given by the 
functions w+(x) and w_(x) respectively, both periodic 
with period L, see Fig.[TJ The particles are subject to an 
external static force f acting along the longitudinal (x) 
direction (potential U(x,y) = — fx). The hydrodynamic 
interactions within the system can be neglected provided 
particles are smalli^ and their concentration is low. We 
define the system of units with m = L = ksT = 1. 
The unit of time in this system is r = L\fmJk^T . The 
velocity v of the particle of particles is then governed by 



with delta-correlated Gaussian random force £ = (£ x , £ y ): 
(&(*)> = 0, (6(^ 3 (s)) = S t3 S(t - s); i,j are x or y. 

The evolution of the probability density function 
(PDF) P (q, v, t) of position q = (x, y) and velocity 
v = (v x , v y ) corresponding to Eq. (fTJ) is governed by the 
Klein-Kramers equation^ dtP = L x>Vx P+ L y Vy P , with 
L q>v = —vd q + d q Ud v + 7o\, [v + d v ] with no-flow condi- 
tion at the boundaries. Approximations to this equation 
can give rise to effective theories concentrating on rele- 
vant x-coordinate and suppressing the irrelevant y-one. 
Let us first discuss necessary conditions for such an ef- 
fective description. 

According to the Bayes theorem, the joint PDF 
of the position and the velocity is given by the 
product P(q,v,i) = $(y,v y \x,v x ,t)p(x,v x ,t) 
of the marginal probability density p(x, v x ,t) = 
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and v v conditioned on x, v x , and t. The fast relaxation 
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dvy P(q, v, t) and the joint PDF of y 
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approximation** assumes that <f>(jj,v y \x,v x ,t) is equal 
to equilibrium PDF of y and v y , conditioned on x: 

$(y,Vy\x,v x ,t) = $(y,v v \x) (2) 

mth*(y,vy\x) = e -l*V*+u{ x ,y)] /^£+(x) e -u( x ,y) dy _ 

In this case the full dynamics, Eq. (JTJ, can be replaced 
by the motion of a particle in the potential A(x) of mean 
force defined by 



dA(x) 
dx 

yielding 



dy J dv y d x U(x,y)$(y,Vy\x) . 

_{x) -oo 



(3) 



1 v x -d x A(x) + y'2 1 i x (t) 



(4) 



The difference between the external force and the mean 
force, SF x (y) = —d x U(x,y) +d x A(x), results in an addi- 
tional effective deterministic force as well as in position- 
dependent effective friction-^ and diffusion coefficients^. 

The assumption Eq. is valid if (i) the distribu- 
tion of y relaxes fast enough to the equilibrium one, 
(ii) equipartition of the kinetic energies corresponding 
to v x and v y holds, and (iii) the two velocity compo- 
nents are uncorrelated at any time. Burada et al.— ana- 
lyzed time scales involved in the problem. These are the 
times t v = 7 Ay 2 /2 and t x = 7 Ax 2 /2 to diffuse over dis- 
tances Ay and Ax, respectively, the characteristic times 
associated with the drift rj rift = min(7 Ax//, 7 Ax/v x ) 
and Tj rift = jAy/vy, and the velocity correlation time 
T cor = 1/7. A general criterion that has to be satis- 
fied is that max (t v /t x , r| r ; ft /Tf rift , r cor /r| rift ) < 1. This 
can be achieved either for strongly anisotropic friction 
Jx 7ij 16 i 18 or for relatively flat boundaries^. 

In the high friction limit, 7 ^> 1, the Klein-Kramers 
equation associated with Eq. (0} simplifies to the FJ 
equation^ dtp(x,t) = d x \e~ A ^d x (e A ( x 'p(x, i))] , and 
the analytic expression for the potential of mean force 
A(x) can be derived along the lines of RefJ^. To do 
this, we mimic the interaction of the particles with 
walls by a quadratic potential growing in the direction 
normal to the wall, U±(n±) = ^n± with interaction 
strength k and n± being the coordinate along the nor- 
mal to the upper or lower boundary taken at the point 
(x,lo±(x)). For x fixed, the energy depends only on y and 
is given by U±(x,y) = 0.5k (y ~ uj±(x)) 2 cos 2 a± with 
a± = arctan (dui±(x) / dx) for y < u;_(x) and y > lu + (x), 
and vanishes otherwise. Eq. (J3j> then reduces to d x A(x) = 
I-00 d y I-oc dv v ( 9 x U + d x U ± ) $(y, v y \x). Integrating 



over v v and interchanging the derivative and integration 



yield d x A(x) 



-d x In 



The 



J^dyexpi-U-U^ 

domain of integration over y can be divided into three 
intervals —00 < y < lu^(x) < y < lu + (x), and 

o; + (a;) < y < 00. Integral J"~ dy . . . can be evaluated 



explicitly by integrating by parts. The corresponding ex- 
pression vanishes in the limit of hard walls, i.e. n — > 00. 
The same happens with the integral J u+ ^ dy . . . . Con- 
sequently, the mean force reads 



dA(x) 



d x In 



j+{x) 



dyexp(~U(x,y)) 



UJ— (x) 



(5) 



The potential of mean force is the free energy 
associated with the partition function Z(x) = 

Iu> + (x) d y ex P ( — U(x, y)) and does not depend on 7. This 
result relies on Maxwell distribution of v y , i.e. on 
equipartition. As we proceed to show, this property 
breaks down if the motion is not overdamped. 

In the following, we study the mobility /i//io = 
lim t _ i . 00 7 (x(t)) /(ft) and the effective diffusion coeffi- 
cient D cfi /D = lim^oo 7 ((x(t) 2 )-(x(t)) 2 )/(2 1) of par- 
ticles moving in a sinusoidally-shaped^ channel with the 
top boundary given by 

uj + (x) = i [Ail + Auj - (AQ - Auj) cos (2tt x)] , (6) 

and with flat bottom boundary oj_(x) = 0,Fig. [T] Af2 
and Auj denote the maximal and the minimal width of 
the channel, respectively, with aspect ratio 5 = Auj/Afl. 

Fig. [5] shows the influence of the external force magni- 
tude / on the particle mobility (see panel (b)) for various 
friction coefficients 7. For 7>1, one observes the known 
dependence of /Vmo on /— . Starting from the asymptotic 
value ^i/^o = 2Vd/(l + 6) for / < 1 (sec Eq. (45) in 
Refi^), the mobility increases with force magnitude till 
the asymptotic value /z/Z-'o = 1 is reached for / — > 00. It 
can be calculated using the Stratonovich formula^ giving 
rise to 



Mo 



P + (27TY 



P + 



(27T) 



(7) 



The result Eq. (JT)) matches perfectly the numerics for 
7 = 100, see Fig. H(b). 

In the diffusion dominated regime, / < 1 in Figs.[2ja)- 
(b), the mobility decreases with decreasing 7. Such a 
dependence is known for arbitrarily damped Brownian 
motion in periodic potentials^ and might witness for 
the applicability of the reduced description. For cho- 
sen geometry [ij [i§ goes as ~ <-y - 5 ; the exponent depends 
both on Af2 and on Auj (not shown). Moreover, the mo- 
bility and the effective diffusion coefficient coincide for 
/< 1 (sec Fig.[2](c)) thus corroborating the Sutherland- 
Einstein relation^. 

For larger forces the mobility increases until it reaches 
its maximum at / ma x (depending on friction) and then 
decays as \ij \i§ <x f~ a with a < 1. The particle mobility 
for the reduced dynamics in the presence of the poten- 
tial of mean force A(x) obtained by simulating Eq. (3]) 
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Figure 2. (Color online) Results of simulation of full dynam- 
ics, Eq. {TJ, in a channel geometry of Fig.[Tjwith Aw = 0.005 
and AQ — 0.05. Panel (a): Particle mobility /u//io as a func- 
tion of 7 for different force magnitudes /. Panels (b) - (d): 
the force dependence of different dynamical characteristics of 
the system for 7 = 1,6 and 100. (b) The particle mobility, 
symbols. Superimposed are numerical results for reduced dy- 
namics, Eq. ([4]) (dashed lines), and the analytical result for 
the overdamped case Eq. (solid line), (c) The effective 
diffusion coefficient, D e s/Do. (d) The 2nd central moment 
a 2 of v v . The horizontal dash-dotted lines indicate unity. 



is shown in Fig. [2] (b) by dashed lines. The approxi- 
mation overestimates the true mobility for all / but is 
sufficiently accurate for 7 > 5 and / < / max - For 7=1 
the discrepancy is large even for /< 1. Introducing the 
position-dependent friction coefficient 7(2:), as proposed 
in Ref^, gives corrections of the order of (Ail) and does 




V - <V > , V - <v > 
x x y y 



Figure 3. (Color online) Normalized probability distribution 
functions of v x and v y as functions of the friction coefficient 
7 for / = 100. The mean values are (v x ) « 3.46,6.78,0.99 
(for 7 = 10°, 10\ 10 2 ) and (v y ) ~ (for all 7 values). The 
black dashed lines indicate the Maxwell velocity distribution 
exp 



{-vl, y /2) /V2^. 



not improve the agreement. 

Our derivation of the effective dynamics implied the 
Maxwell distribution of v y and the homogeneous distri- 
bution of y. In Fig. [3] we present the velocity distri- 
butions P(v x ) and P(v y ) centered at their means (v x ) 
and (v y ) for different 7 and for fixed / = 100. The 
distribution P(v x ), Fig. [3] (a), undergoes a transition 
from a normal (Maxwell) distribution with variance 1 for 
7 = 100 over a broader bimodal distribution for 7 = 10 
to the broad symmetric function for 7=1. In con- 
trast, P(v v ), Fig.[3](b), stays symmetric regardless of the 
value of 7. Similar to P{v x ), the smaller the friction the 
broader the distribution. The deviation of P(v y ) from 
the Maxwell distribution is given by its second central 
y ) 2 , Fig.[5](d), which is unity (in 



moment o" 2 = (vfy 
our units) for the Maxwell one. The distribution of v y is 
Maxwcllian independent of / for sufficiently high friction 
for any /. For smaller friction the transversal velocity 
distribution broadens, a 2 cx f@ with j3 > 1, if / exceeds 
a critical magnitude f c which practically coincides with 
/max- In other words, the decrease of the particle mo- 
bility goes hand in hand with violation of equipartition. 
The 2nd central moment of v x shows the same depen- 
dence on the external force magnitude and on the friction 
coefficient as the one of v v . 
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Figure 4. (Color online) Particle trajectories for 7 = 1 (black 
solid line and dashed line) and for 7 = 100 (bright solid line). 
The force magnitude is / = 100. The geometry parameters 
are AQ = 1 and Auj — 0.1. 



Finally, we discuss the impact of 7 on the effective 
diffusion coefficient. In the high friction limit we re- 
produce the results of Rcf^: starting from the value of 
D cS /D = 2VS/(1 + S) for / ->■ the effective diffusiv- 
ity grows with increasing / until it reaches its maximum 
and then decays and finally tends to D c ^/D a = 1 for 
/ — > 00. For / -c 1 the effective diffusivity grows with 
friction. For sufficiently strong forces, / > 10, the be- 
havior of D c ff/Do as a function of 7 is non-monotonic, 
passing through a maximum at some value of 7 which 
depends on the force magnitude and on the channel's 
geometry. 

Let us discuss the nature of the equipartition viola- 
tion. In free motion the velocities "forget" their initial 
values and assume equilibrium distributions for t > I/7 
for any 7. The value of 7 influences however the proper- 
ties of confined motion, see Fig. [4J changing from erratic 
for high friction to almost regular - known for the deter- 
ministic case^S - for 7 —> 0. Due to the reflection at the 
boundaries the acceleration caused by the external force 
/ is transferred from the longitudinal to the transversal 
velocity component violating equipartition and leading to 
the monotonous growth of a 2 with /. The behavior in the 
x-direction is more complex. On one hand, the particles 
reflected at an "optimal" angle can fly over several cells 
to the left or to the right. We found that the probabil- 
ity of long stretches increases with force magnitude and 
with decreasing friction, and that long excursions into 
the direction of the force are always more probable than 
in the opposite one. On the other hand, some particles 
undergo repeated collisions with walls where these are 
almost parallel, i.e. close to the minima and maxima of 
the channel's width, see dashed line in Fig. [4] leading to 
trapping. The interplay of both effects leads to complex 
behavior and non-monotonicities described above. 

Conclusions Let us summarize our findings. We in- 
vestigated the impact of the friction coefficient on bi- 



ased Brownian motion of point-like particles in a two- 
dimensional channel with smoothly varying width. For 
high friction, the adiabatic elimination of the transverse 
degrees of freedom results to an effective description for 
the slow ^-coordinate involving the potential of mean 
force A(x), leading to an exact analytical result for 
the particle mobility Eq. |[7J) valid for narrow channels^. 
Comparing the results for reduced dynamics with numer- 
ical results for the full problem, we show that the reduced 
description overestimates the mobility, although the ac- 
curacy of the approximation is sufficiently good for 7 > 5 
and small forces / < / max . There exists a characteristic 
force strength f c beyond which the reduced description 
fails. The force f c become less with decreasing friction. 
The origin of the failure of the effective description is the 
violation of equipartition for the fast coordinate y and 
velocity v y . 

The authors arc thankful to P. Hanggi and G. Schmid 
for useful discussions and acknowledge financial support 
by the VW Foundation via project 1/83903. 
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